Maximal Slicing for Puncture Evolutions of 
Schwarzschild and Reissner-Nordstrom Black Holes 
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We prove by explicit construction that there exists a maximal slicing of the Schwarzschild space- 
time such that the lapse has zero gradient at the puncture. This boundary condition has been 
observed to hold in numerical evolutions, but in the past it was not clear whether the numerically 
obtained maximal slices exist analytically. We show that our analytical result agrees with numerical 
simulation. Given the analytical form for the lapse, we can derive that at late times the value of the 
lapse at the event horizon approaches the value jg\/3 w 0.3248, justifying the numerical estimate of 
0.3 that has been used for black hole excision in numerical simulations. We present our results for 
the non-extremal Reissner-Nordstrom metric, generalizing previous constructions of maximal slices. 
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I. INTRODUCTION 

When decomposing the 4-dimensional Einstein equa- 
tions into (3+l)-dimcnsional form, a crucial role is 
played by the lapse function that determines how the 
4-dimcnsional manifold is sliced into 3-dimcnsional hy- 
persurfaces. For spatial hypersurfaces the lapse deter- 
mines how time progresses between hypersurfaces. There 
is large freedom in choosing a lapse, but a particularly 
interesting possibility is to determine the lapse from the 
condition that the mean extrinsic curvature of the slices 
vanishes at all times. This condition corresponds to a 
certain maximal volume property of the hypcrsurface [1] , 
and the resulting gauge condition is referred to as maxi- 
mal slicing. 

Maximal slicing has played an important role in both 
analytical and numerical work using 3+1 decomposi- 
tions. It was suggested as geometrically motivated anti- 
focussing slicing condition by A. Lichncrowicz [2] al- 
ready in 1944, i.e. maximal slices avoid certain coordinate 
pathologies. Furthermore, maximal slicing helps to deal 
with physical singularities associated with black holes. 
For a single Schwarzschild black hole, maximal slicing 
can be constructed analytically [3, 4]. One can obtain 
maximal slices that extend from the "right-hand" spatial 
infinity to the "left-hand" inner infinity of the extended 
Schwarzschild spacetime, and these slices give a complete 
foliation outside the event horizon while approaching a 
limiting slice inside that does not reach the physical sin- 
gularity. In this sense maximal slices avoid the physical 
singularity, and even in more general situations maximal 
slices have been found to be singularity avoiding [5]. 

For these reasons, maximal slicing has been used fre- 
quently in numerical simulations of one black hole, e.g. 
[3, 6, 7], and also in binary black hole mergers in full 3D, 
e.g. [8, 9]. The lapse for maximal slicing is obtained in 
general by solving an elliptic equation, but even though 
there now are algebraic lapse conditions that are com- 



putationally far less expensive and share the singularity 
avoidance property, maximal slicing is still being used 
when smoothness of the lapse is an issue, say in the 
currently most advanced computations of gravitational 
waves from black hole mergers [10]. In addition, there 
are shift conditions [11] which in many cases overcome 
the so-called slice stretching problem of maximal slicing 
that previously limited black hole evolutions like [9] to 
short evolution times. 

In this paper we establish the connection between the 
analytically known solutions for maximal slicing of the 
Schwarzschild spacetime [3, 4] and the maximal slicing 
computed in the puncture evolution method for black 
holes [8]. In the puncture 1 method, the left-hand asymp- 
totically flat region of a black hole is effectively compact- 
ified by analytically factoring out a coordinate singular- 
ity and working on R 3 , which can result in a significant 
technical simplification over approaches which implement 
an inner boundary for black hole excision. However, al- 
though the general solution to the linear elliptic equation 
for the maximal slicing lapse is known for Schwarzschild, 
this solution depends on the boundary condition imposed 
on the lapse. Only "odd" and "even" boundary con- 
ditions have been considered before, and the numerical 
lapse of the puncture evolutions is not of this type. Fur- 
thermore, there was no rigorous proof that in the pres- 
ence of the puncture coordinate singularity a regular so- 
lution to the lapse equation exists, and while numerically 
quite robust, it was not clear that the numerics was able 
to reliably determine this lapse. 

The key idea in this paper is to impose the numerically 
observed behaviour of the lapse, namely that its first 
derivative vanishes at the puncture, as a novel bound- 
ary condition on the lapse. This boundary condition we 
call the "zero gradient at the puncture" ("zgp") condi- 
tion. Since it is the coordinate singularity at the puncture 
that implies the "zgp" condition, we refer to the resulting 
maximal slicing lapse also as the puncture lapse. 

The main result of the paper is that we can prove by 
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explicit construction that there exists a maximal slic- 
ing of the Schwarzschild spacetime such that the lapse 
has zero gradient at the puncture, and furthermore that 
our analytical result agrees with numerical simulation. 
Concretely, we obtain the "zgp" or puncture lapse as a 
time-dependent linear superposition of the odd lapse and 
the even lapse. As an application of the analytical form 
for the puncture lapse, we can derive that at late times 
the value of the lapse at the right-hand event horizon 
approaches the value 0.3248, justifying the nu- 

merical estimate of 0.3 that has been used for black hole 
excision in numerical simulations. 

There are two technical difficulties to overcome. First, 
the analytic formulas of [3, 4] are non-trivial, in particu- 
lar the calculation is based on Schwarzschild coordinates 
which lead to complications in the numerical evaluation 
of the resulting integrals. And second, we have to estab- 
lish an explicit mapping between the coordinates used 
in our analytical study and the numerical simulation be- 
fore we can perform a comparison. In this paper we only 
consider the case of vanishing shift. 

Since we have to discuss [3, 4] in quite some detail, we 
take the opportunity to derive all relevant equations not 
only for Schwarzschild, but for a non-extremal Reissner- 
Nordstrom black hole. This turns out to be possible with- 
out major changes to the formalism and extends in some 
places previous results [12]. 

The paper is organized as follows. In Sec. II, we recall 
the analytic solution to the lapse equation. In Sees. Ill 
and IV, we discuss in some detail the odd and even lapse, 
respectively, since these solutions will be superimposed to 
obtain the "zgp" lapse. Furthermore, in Sec. IV we make 
explicit contact between the equations given in [3] and 
[4]. In Sec. V, we construct the "zgp" lapse and discuss 
the late time limit of the lapse at the right-hand event 
horizon. In Sec. VI, we compare the analytical results 
for odd, even, and "zgp" maximal slices. In Sec. VII, we 
construct the coordinate transformation needed to make 
contact with the numerical "zgp" lapse and compare with 
our analytical result. We conclude in Sec. VIII. 



II. LAPSE FOR MAXIMAL SLICING OF THE 
REISSNER-NORDSTROM SPACETIME 

To derive a differential equation for the lapse which 
maximally slices the Reissner-Nordstrom spacetime, we 
follow the reference [3] closely. There in 1973 - concen- 
trating on even boundary conditions - F. Estabrook and 
H. Wahlquist were the first to find and S. Christcnscn, 
B. DeWitt, L. Smarr and E. Tsiang to verify numer- 
ically the analytic solution for maximal slicing of a 
Schwarzschild black hole. However, to simplify some cal- 
culations for a comparison of different boundary condi- 
tions imposed on this lapse, as done in Sec. Ill to V, 
the notation of [4] will be used to some extent and will 
be generalized to include electric charge. In the latter 
paper, R. Beig and N. O Murchadha rc-derived this fo- 



liation based on a more geometrical ansatz and studied 
for even boundary conditions the late time behaviour of 
the throat in the limit as proper time at infinity becomes 
arbitrarily large. We will make contact between formulas 
of those two references in Sec. IV. 

We start with the source-free Einstein-Maxwell equa- 
tions 

G>„ := i? M „ — — Rg^ = SnTftv, (1) 
VF^ = 0, (2) 

for a charged black hole, where the stress-energy tensor 
Tfiu is given in terms of the Maxwell tensor by 

Tfu, — -^{FnaFJ 7 — ~79tu>F ae F ae ). (3) 

The 3+1 split of the equations in the ADM-formalism 
(see R. Arnowitt, S. Dcser and CW. Misner [13] or [1]) 
yields the Hamiltonian and momentum constraints 

H:= R + K 2 - KijK ij - 16irp = 0, (4) 
M z := V J (^- ? '-7^A')+8 7 r Sl =0, (5) 

together with the evolution equations 

(dt-Cp)'fij = -2aKij, (6) 
{dt-CffiKy = {-V i V j + R ij +KK ij -2K ik K k j )a 
-(8irSij - 4n(S - phij) a. (7) 

Here the fundamental dynamical variables are the in- 
duced, positive definite 3-metric 7y and the extrinsic 
curvature Kij. Note that Cp is the Lie derivative with 
respect to the shift vector and V* is the covariant spa- 
tial derivative associated with jij. Furthermore, Rij is 
the 3-dimensional Ricci tensor and the Ricci scalar R is 
its trace. The stress tensor SV;, the momentum density 
vector Si, and the total energy density p are obtained as 
projections of the Maxwell stress-energy tensor along 
the normal n M by 

p = T^n^rf , Si — T v 7i M , and .S',, /,,. (8) 

The mean extrinsic curvature 

K := jijKV = -V^ (9) 

measures the amount of "crunch" , i.e. the fractional rate 
of contraction of 3-volume along a unit normal to the 
surface. 

In the following we want to restrict ourselves to maxi- 
mal slicing, 

K = 0. (10) 

The reason for its name can be inferred from a variational 
principle maximizing the volume 

V(n) = f Vdet { lmn } d 3 x (11) 
J n 
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of a bounded but arbitrary portion 51 of the Cauchy 
slice £, where as pointed out by J.W. York in [1] it follows 
that the trace of the extrinsic curvature has to vanish. 
Contracting the evolution equation for Kij, (7), in the 
context of K = one obtains the maximality condition 



Aq := V'Vid = (R + 4tt(5 - 3p)) a, 



(12) 



which is an elliptic equation for the lapse function. 

As in [3, 4, 12] the calculations are performed in the 
radial gauge 
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respectively. Furthermore, the ##-componcnt of the evo- 
lution equation of the extrinsic curvature, (7), is used 
together with the maximality condition (12), which for 
the radial gauge reduces to a condition involving the rr- 
and ##-componcnt of (7) only, to obtain as further equa- 
tions 
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ds 2 = {-a 2 + — ) dr 2 + 28 drdr + 7 dr 2 + r 2 dtf (13) 

7 

with a, 3 and 7 being functions of r and r only. Whereas 
{t, r, 9, <p} are the standard Schwarzschild coordinates 
with the Schwarzschild radius r measuring the circum- 
ference divided by 2ir or the square root of area divided 
by Air of the 2-sphere, the function r = r(t, r) is clearly 
not. 

Note furthermore that for the Maxwell tensor F^ v of a 
purely radial electrostatic field one can make the ansatz 
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(20) 



0. (21) 



The equations (18) to (21) generalize the Schwarzschild 
black hole present in Eq. (2') to (5') of [3] to a black 
hole carrying electric charge. (Note the correction of two 
signs in [3].) 

The momentum constraint (19) immediately yields 



see e.g. Chap. 18 of [14], and one can readily verify that 



En = 



(15) 



solves the source- free Maxwell's equation (2). For (static) 
Schwarzschild coordinates the classical result Eq = ^ 
for the electric field of a point particle with charge Q 
situated at the origin is recovered. Finally, the stress- 
energy tensor can be calculated according to (3) and its 
projections (8) are 



/3(r,r) 



u(t, r)7(r, r) 



C(r), 



(22) 



where C is a function of r only. This equation can be 
used to replace 3 in the remaining equations such that 
from the Hamiltonian constraint a first-order partial dif- 
ferential equation (PDE) for 7 arises. Its solution is given 
by 



7(7-, r) = 



1 



(23) 



P = 



87rr 4 ' 



= 0, and Sij = p diag [■ 



-7,r 2 ,r 2 sin 2 9] 



In addition, we observe that the trace of SV 



(16) 
coincides 



with p, S := r ) lJ = p. 

One can readily show that the equation for the vanish- 
ing of the trace of the extrinsic curvature, K = 0, can as 
in Eq. (1') of [3] be written in the form 



(log7) >T = ^ ( 



/3V 



7 



(17) 



which can be used in the following to eliminate deriva- 
tives of 7 with respect to t. After some rewriting, the 
Hamiltonian and momentum constraints, (4) and (5), 
yield 



.7 



a 2 7 



■7(1 



1 

1 



Q' 2 



(18) 



where using the function 

r H 
for convenience the polynomial 

p c (T,r) = r 4 f(r) + C 2 (T) 

= r 4 -2Mr 3 + Q 2 r 2 + C 2 (r) (25) 

has been introduced. Note that M, which arises as a 
constant of integration, could in principle be a function of 
r, however as in [3] it follows from the evolution equation 
(20) together with (17), (22) and (23) that M must be 
r- independent. It is worth pointing out that the radial 
metric function we obtain in (23) generalizes Eq. (9) of 
[3] by the charge, and the same 3-metric for maximal 
slicing of Rcissncr-Nordstrom as found by M. Duncan in 
[121 is obtained. 
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Finally, making use of j3, (22), and 7, (23), the maxi- 
mal slicing condition (17) reduces to 



f^d£ 
r 2 dr ' 



(26) 



which is the differential equation for the lapse we were 
looking for. By integration its general solution is found 
to be 



a(r, r) 



% / 

dr J r 



V A dy 
Pc(r,y)i 



(27) 

where C and D are functions of t only and one inte- 
gration limit has been fixed without loss of generality 
at spatial infinity. Three different boundary conditions 
leading to an odd, even, and puncture lapse will be dis- 
cussed next. 



III. ODD LAPSE 
A. Derivation 

The simplest special case of the lapse (27) is, because of 
its underlying antisymmetry with respect to the throat, 
referred to as odd (in some references also as antisym- 
metric) and corresponds to the boundary conditions 



lim 



odd 



= ±1 Vr G , 



(/(/■ 



(28) 



To measure proper time at infinity (on the right-hand 
side of the throat denoted by a superscript "+" ) the lapse 
is unity there and implied by the antisymmetry the lapse 
is minus one at the puncture (on the left-hand side and 
hence denoted by a "— "). The odd lapse 



odd 



(C,r) 



dC 



dr,, 



(29) 



0. The 



is obtained by setting D{T ot id) = ±1 and 
time-independent slice label C can be chosen indepen- 
dent of the time at infinity r a dd and is, as pointed out 
for Schwarzschild in Appendix A of [4], purely gauge for 
Reissner-Nordstrom. 



B. Schwarzschild lapse 

The outer event horizon and inner Cauchy horizon at 
r± = M ± y/ M 2 — Q 2 are given by the real roots of /(r), 
which are time independent, while the throat rc is ob- 
tained as a real root of pc (r) and therefore depends on 
the choice of C. Hence, starting with the throat in the 
same location as the event horizon implies C = 0. In this 
case one obtains from (29) the so-called Schwarzschild 
lapse 



odd 



(C = 0,r) = ±v^W, 



(30) 



from (23) the radial component of the metric 
7(r) = 



f(r) 



(31) 



and from (22) a vanishing shift. Hence the static 
Reissner-Nordstrom metric written in Schwarzschild co- 
ordinates, 



ds 2 = -f(r) dt 2 + 



f(r) 



dr 2 



'*dVL 2 



(32) 



is recovered identifying Q with the charge and M with 
the mass of the black hole (0 < \Q\ < M). 

This result was of course expected, since it is well- 
known that with this coordinate choice the spacelike 
slices given by t = const, being perpendicular to the time- 
like Killing vector in the exterior regions, arc maximal. 
That one can immediately sec from the definition of the 
extrinsic curvature Kij, (6), which for zero shift reduces 
to 



(33) 



Since the 3- metric 7^ is time-independent, each of the 
components and therefore also the trace of the extrin- 
sic curvature vanishes, hence K = 0. Note furthermore 
that by contraction of (33) setting K = one immedi- 
ately obtains the statement that for maximal slicing with 
vanishing shift the determinant of the 3-metric is time- 
independent. Hence the singularity avoiding property 
comes to light as the variation of the local volume re- 
mains fixed [15]. 



C. Isotropic lapse 

By applying for vanishing shift a purely spatial coor- 
dinate transformation r = r(R), the hypersurfaces obvi- 
ously remain maximal. To obtain isotropic coordinates 
by writing dr 2 = {-j^) 2 dR 2 and by comparing the dt 2 -, 
dr 2 - and d51 2 -terms the set of equations 



a&|(C = 0,J2) = f{r{R)), 



* 4 (i?) 



1 



( — ? 



f(r(R)) y dR 
R 2 ^(R) = r(R) 2 , 



(34) 
(35) 
(36) 



is read off. From (35) and (36) the ordinary differential 
equation (ODE) 



dr 



r(R) 



dR = ± ^m R 



(37) 



can be inferred where r can take values r > r + only. By 
integration the solution is found to be given by 



r(R) 



2R' AR 2 



R + M 



M 2 - Q 2 
4~R 



(38) 
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choosing an arbitrary constant of integration such that 
r and R coincide at infinity. With R + = ^^/M 2 — Q 2 
corresponding to r + this relation can also be piecewise 
inverted to yield 



R(r) = f |(-M + r(l-V7W)) 
n[r) 1 if-M + rfl + yjjtr))) 



for 



< i? < R+ 
R > R+. 
(39) 

From equation (34) follows the isotropic lapse 



„± rr-n m (^r + m)(2R - m) + q 2 
<* odd (c - o, R) - (2i? + M)2 „ Q2 (40) 

and from (36) one can read off the conformal factor 

M , n 



^(R) = (1 



2R' 



4i? 2 J 



(41) 



Hence the Rcissncr- Nordstrom metric in isotropic coordi- 
nates as found e.g. in [16] is recovered which for vanishing 
charge reduces to the well-known Schwarzschild case. 



D. Isometry 

For the case < \Q\ < M, one can see from the co- 
ordinate transformation (38) that there is an isometry 
present, since the value of r and therefore also the line 
element (32) remains the same under a mapping 



R 



M 2 - Q 2 
AR 



(42) 



These two values of R for every value of r > r + result 
in two isometric parts of the spacetime (in terms of 
Kruskal-Szekeres diagrams the regions I and I' in Fig. 3 
for Schwarzschild). In particular, R = is simply a com- 
pactificd image of infinity in the other universe referred 
to as puncture. The fixed point set at the event horizon 
r + is a minimal 2-sphere of an Einstein-Rosen bridge [17] 
located at R + , the so-called throat (see the embedding 
diagram of Fig. 2). 

Note also that the odd lapse in isotropic coordinates, 
(40), changes its sign with positive values on the right- 
hand side of the throat being the original unextended 
space and negative ones on the left-hand additional re- 
gion obtained by analytic extension of the spacetime. 
In particular, the odd lapse vanishes on the throat, i.e. 
a ± (C = 0, i?+) = (c.f. the lapse profile as shown for 
Schwarzschild in the upper plot of Fig. 6). 



IV. EVEN LAPSE 

A. Derivation I (via Einstein-Maxwell equations) 

The even lapse is symmetric with respect to the throat 
and its boundary conditions are usually stated as unit 



lapse at infinity and a vanishing gradient at the throat. 
Since by symmetry time has to run equally fast at both 
spatial ends of the manifold, one can also formulate even 
boundary conditions as 



lim 



(43) 



This motivates for the differential equation obtained for 
the lapse in (27) the particular choice D(r even ) = 1 made 
by F. Estabrook et al. and generalized here, so 



i 



dC 



V 4 dy 



\ dr even J r pc{yY' 

(44) 

Since one starts with the 3-geometry of the t = T even = 
hypersurfacc, i.e. with the initial radial component of 
the 3-metric given by 7 = jct, one can infer from (23) 

that C{T eV en 

= 0) = initially. The function C(T even ) 
though is still undetermined as only single sheets of 
the 3-gcomctry labeled by C have been looked at so 
far. As stated in [3], the r et , erl -dependcncc of C can be 
fixed by imposing the requirement of smoothness across 
the Einstein-Rosen bridge by passing to the Reissner- 
Nordstrom line element (32) in Schwarzschild coordi- 
nates. For the "height function" t = t(C(r even ), r), com- 
paring dr 2 ven -, dT even dr- and c?r 2 -terms of the metric the 
following two PDEs, 



m 



dr e . 



at 



C 



9r /Mv^W : 



are found. Making for (46) the ansatz 



i(C,r) = 



Cdy 



f(y)Vpc{y) 7 



(45) 
(46) 

(47) 



one finds that this equation also satisfies (44) and (45) for 
sc coinciding with rc- The latter is the unique (double 
counting for C = Cu m ) larger real root of the polynomial 
Pci r ) m [ r Cu m > r +] an d determines the radial coordinate 
of the throat on a given maximal slice. Hence 



teven(C,r) = H C {r), 

where we have introduced the integral 



H c (r) := - 



Cdy 



f{y)y / Pc{y) 



(48) 



(49) 



for r >rc- Note that in (49) the integration across the 
pole at r+ is taken in the sense of the principal value, 
and by the same arguments as in [4] the corresponding 
slices extend smoothly through the event horizon r + and 
throat rc- 

The procedure of deriving (48) with (49) is non- 
trivial as pointed out in detail in [3]: The require- 
ment of smoothness across the bridge is t even — > and 
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teven —> oo as r approaches its smallest value rc at 
the center of the bridge, or in terms of Kruskal-Szekeres 
coordinates -^T = as X — > 0. As mentioned in [18], 
this fixes the value of C(rc) such that pc{ r c) = holds. 
In particular, with rc becoming at infinite times a double 
counting root rc Um with the value 



1 



ro Um = 7 (3M + V9M2 



J ) >0 



corresponding to 



C; 



lim 



oo) 



(50) 



(51) 



C(r evt 



the singularity avoidance of the maximal slicing of the ex- 
tended Reissner- Nordstrom spacetime as proven in corol- 
lary 3.8 in [5] becomes clear. 

Finally, since t even (r) and r even have to coincide at 
spatial infinity modulo an arbitrary constant, which can 
be set to zero, from (48) in the limit r — > oo one can infer 



Teven(C) = lim t even (C,r) = He (oo). (52) 



We want to stress that in (47) the particular choice of 
the larger root of pc(f) has been made since we are in- 
terested in the "horizon-horizon" subfamily of maximal 
slices only, which extend to both spatial infinities. Here 
the slices start with the initial hypersurface representing 
the 3-geometry at the moment of time symmetry when 
the throat is at the event horizon r + and develop in- 
ward to approach rc Hm asymptotically and hence avoid 
the singularity. As pointed out in remark 8 of [3] and 
in Sec. Ill of [12], a second "singularity-singularity" sub- 
family of maximal slices can be obtained choosing the 
unique (double counting for C = Cu m ) minor real root 
fc of pc (r) in [r_,rc lim ] instead. The initial hypersur- 
face is now the time-symmetric slice with the throat at 
the Cauchy horizon r_ and the slices intersecting the 
singularity at r = grow outward to the limiting sur- 
face rc Hm ■ As discussed in the next subsection, one can 
compare the situation with the mechanical analogue of a 
particle of energy E = C 2 travelling in a repulsive poten- 
tial V(r) = — r 4 /(r), where rc and fc represent turning 
points (c.f. Fig. 1). Thinking of particles with energies 
ranging from E = to E lim = V(rc Um ) = Cf im it is clear 
that both subfamilies together cover the spacetime with 
spacelike maximal slices (c.f. the Kruskal-Szekeres dia- 
gram for Schwarzschild and the Carter-Penrose diagram 
for Reissner-Nordstrom in Fig. 1 of [3] and [12], respec- 
tively). For its limited use in numerical relativity, how- 
ever - as the slices lie always completely inside the event 
horizon and hit the singularity - we will not investigate 
the "singularity-singularity" subfamily further. 



B. Contact between [3] and [4] by a mechanical 
analogue 

To make contact between formulas found by 
F. Estabrook et al. [3] and those derived by R. Beig and 
N. O Murchadha [4], one can apply an energy conserva- 
tion equation worked out in Appendix B of the latter 
paper. There two functions F and J are introduced by 



F(E,r) = f 

J XT 



W(y) dy 



(53) 



J(E,r) = / ^E - V(y)V(y)W(y) dy. (54) 



In the following, W may have a simple pole at x = x 
and V mapping xe to E and a; to is a smooth function 
V : [x , oo [ — ► R satisfying < E < V(xq), V(x) < E, 
and V'(x) < for x > xq. By differentiating J twice with 
respect to E the formula 



1^1 

dE 2 



J{E,x) = - 



V(x)W(x) 



VE - V{x) V'{x) 



y/E - V(y) dy 



V(y)W(y) 



V'(y) 



F(E,x)+2E—F(E,x) 



dy 
(55) 



is derived, c.f. Eqs. (B7) and (B9) in [4]. 

Now, as pointed out in [19] using proper velocity 4^, 
i.e. the rate of change of proper time along a maximal 
hypersurface with proper distance, from 



[ J 



r 4 f(r)=T + V (56) 



an energy conservation equation arises for a particle of 
total energy E = C 2 moving in the repulsive potential 

V = -r 4 f(r) with a kinetic energy T = C 2 (j^f . By 
demanding 4^ = for unbounded particles starting from 
infinity the closest approach to the singularity - or for 
bounded ones starting from the singularity the outermost 
radial distance - is found to be given by the two real 
roots rc and fc of the polynomial pc{r), c.f. Fig. 1. 
Furthermore, R. Beig and N. O Murchadha [4] point out 



that for E > Eu m with En 



V(r Glim ) = Cf im maximal 



slices hitting the singularity at r = are obtained [53]. 
They also argue that F(E,r) can be thought of as the 
time it takes a particle of energy E to travel from rc 
to r. In particular, in an analysis following in a later 
paper [20] based on [4, 21] we will be interested in the 
way this function blows up as E approaches with Cf im 
the maximum of the potential in order to study the late 
time behaviour of the slices. 
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FIG. 1: For maximal slicing of Schwarzschild the mechanical 
analogue consisting of a particle of energy E = C 2 travelling 
in a repulsive potential V = —r 4 f(r) is shown. For the energy 

f "c 2 o "I 
levels E — < 0, — s 52 *-, Cf im > the innermost and outermost ra- 
dial distances for unbounded and bounded particles denoted 
by °rc and fc, respectively, can be inferred. 

One can readily check that the identifications 

W(r) = -jfa, V(r) = -r 4 /M, E = C 2 (57) 

and 

r o = rc Hm , r E = r c , r = r+ (58) 

satisfy the conditions on W, V, and E as stated above 
also in the Rcissncr-Nordstrom case and lead to 



H c (r) = CF(C 2 ,r). 



(59) 



Hence for the derivative of Hc{r) with respect to C the 
important formula 

r 2 1 

= 7T2 ; -K c (r) (60) 

2(r-|M+f) v ^) 2 

can be derived by straightforward calculation using (55), 
where for r > rc we have introduced the integral 



K c (r) 



y(y-3M+^-)dy 



(y-lM + ^) 2 Vp^y) 



(61) 



Note that Kc{r) diverges in the limit C — ► Cu m as rc Hm 
becomes a double root of pc( r ), the square root of which 
appears in the denominator of the integrand. Note fur- 
thermore that from (60) in the limit r —> oo the equation 



—H c (oo) = 
is found, c.f. Eq. (B12) in [4]. 



\ K c (oo) 



(62) 



In order to write the expression (44) for the even lapse 
in terms of the integral Kq (r) , one can readily verify that 

^(U(C,r)-r e „(C)) = — (H c (r) - H c (oo)) 



oo 4 



v d y 

holds and hence the even lapse can be written as 



a toen( C > r ) 



Vp^r) dC dt 
r 2 dr even dC 

if 1 



(63) 



(64) 



K c (oo) yr-^M+^r 



which generalizes Eq. (2.23) and (2.24) in [4] by 
the charge. By the same arguments as applied for 
Schwarzschild in this reference, the term d ^ c dt ^ n in 
the first line of equation (64) blows up at the throat, i.e. 
for r — ► rc, in such a way that the lapse has a smooth 
limit there. To our knowledge the lapse function for- 
mulated in (44) or (64) for maximal slicing Rcissncr- 
Nordstrom (i.e. with electric charge included) has not 
appeared in the literature yet. 

For Schwarzschild, the even lapse profile in terms of 
the Schwarzschild radial coordinate r for the times at 
infinity r even = {0, IM, 2M, AM, 8M } is shown in Fig. 6. 

Recalling that the shift is given by (22) and the radial 
component of the metric by (23), with 



fteven (C? 



Qeven (Cj 7*)^feven (Cj f) 



c 



and 



7euen(C, r) 



f(r) + Pc(r) 



(65) 



(66) 



the full 4-mctric is found. 



C. Derivation II (by a variational principle and/or 
via construction from Reissner-Nordstrom metric) 

In this Subsection we present two alternative deriva- 
tions of the maximal slices, which the reader may skip 
without loss of continuity. 

Note that the maximal slicing condition K = fixes 
the PDE for J£ as in (46) only, whereas boundary condi- 
tions have to be specified to obtain t(C, r) by integration. 
The latter can always be written as the sum of the inte- 
gral Hc(r), (49), and a time translation function of C, 
fixed initially by demanding that for C = one starts 
with the time-symmetric t = hypcrsurface, and deter- 
mined during the evolution by boundary conditions. In 
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particular, for even boundary conditions the throat rc 
has to stay at t = because of the symmetry, hence with 
Hc{ r c) — it follows that this additional function of C 
is identical zero and (48) is found. A non-trivial example, 
namely "zgp" boundary conditions leading to t(C, r) as 
in (97), will be discussed in Sec. V. 

To emphasize the importance of the PDE (46) gov- 
erning the foliation of Reissner-Nordstrom into maximal 
slices formulated as level sets of 



a = t — t(C, r) = const, 



(67) 



where C is a parameter depending only on the new time 
coordinate t, two further ways of re-deriving this equa- 
tion will be outlined here. 

M. Duncan based on a variational principle shows 
in [12] that for Reissner-Nordstrom maximizing the 
3- volume of the line clement 



^ds 2 



-f(r) 



1 



dr 2 



r 2 dfl 2 (68) 



the Euler-Lagrange equation 



D 



' — 



Dr 



i = 



DY 



dY 
dX 



dt dY 
dX dt 



dY 
' dt, x 



(69) 



is ob- 



using the notation OT - ^ , , ^ 

tained. Eq. (69) is satisfied if the expression within the 
square brackets is a function of r only and one may read- 
ily check that (46) is recovered by setting this function 
to C(r). 

Another way followed for Schwarzschild in [4, 22] for 
maximal and in [19] for constant mean curvature slices is 
to examine the behaviour of the normal 



n t dt + n r dr = TVV (t - t(C, r)) = N [ dt 



dt 
dr 



dr 



(70) 

where with the underlying Reissner-Nordstrom line ele- 
ment (32) the normalization constant TV is fixed by de- 
manding 




= TV 

Note that TV is also given by 

TV 



= -1. 



" 1 dt 



(71) 



(72) 



where is the static Killing vector. Hence N can 
be viewed as the boost function [22] of Jj relative to 
a = const. With (9) it follows that the trace of the ex- 
trinsic curvature can be written as 



K = -V M n" = 



1 d [r 2 n r \ 
r 2 dr 



dr 



(73) 



which for maximal slices, K = 0, yields again (46) by the 
same argument as for (69). 

One may readily verify that from the normalization 
(71) it follows that TV coincides with the odd lapse, so 



TV = ± 



"odd- 



Furthermore, writing ^ as 



dt 



= TVn M + with = 



(74) 



(75) 



and using (72) together with n p 



TV#. 
dr 



-aV(,T, one can find 



(76) 



In the context of even boundary conditions this yields 
as in the first line of (64). Here a dd 



Oiodd 



dt cve 



and a eV en are linearly independent solutions of (12) as 
pointed out in detail in [4, 22]. 

One should furthermore mention that instead of vary- 
ing the parameter C with time - which is the choice lead- 
ing to a foliation, i.e. a "time evolution" , made in the lat- 
ter references and in our paper - one could also as done 
in [19] consider adding constant "time translations" to 
t(C, r) for fixed C. 



D. Initial lapse profile 

In order to calculate the initial lapse profile note that 
for C = the primitive for Kc(r) as in (61) is given an- 
alytically by 



K c =o(r) 



2r{2Mr-Q 2 )^J(r) 



{M 2 - Q 2 )(2r 2 - 3Mr + Q 2 )' 
and in particular in the limit r — ► oo by 

2M 



-ftTc=o(oo) = lim K c =o{r) = - 



M 2 -Q 2 ' 



(77) 



(78) 



With (64) therefore the initial lapse profile in 
Schwarzschild coordinates is found to be 



l (C = 0,r) = l 



Mr 



(79) 



taken symmetrically across the throat with the value 
[ (C = 0,r + ) = y / l 



M- 



there. In terms of the 



Schwarzschild isotropic coordinates as in (38) with the 
isometry (42) present this can be written as 



XC = 0,R) = 1- 



AR 



M 4R 2 + AMR + M 2 -Q 2 



(80) 

Obviously, in the Schwarzschild case azf ven (C = 0,r) = 1 
is a solution of (12), while with increasing charge the 
source term 47r(S' — 3p) in the maximality condition leads 
to an increasingly collapsed lapse profile. 
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V. PUNCTURE LAPSE 

A. Motivation of the zero gradient at the puncture 
("zgp") boundary condition 

As stated in the introduction, the motivation to define 
and study the "zgp" boundary condition derives from 
puncture evolutions of black holes. In any numerical 
black hole evolution one has to specify how the black hole 
singularity is treated. In particular, one can excise the 
black hole from the domain at or inside the event horizon, 
thereby creating an inner boundary of the numerical grid, 
one can use the throat as an inner boundary and impose 
isometry conditions there, and finally one can work on a 
punctured R 3 , i.e. R 3 minus one point (the "puncture") 
for every asymptotically flat end inside a black hole. 

Under different names, punctures were studied as early 
as [23, 24] for Schwarzschild and axisymmetric black hole 
data. The puncture topology for black hole initial data 
was revived in [25] as a simple method to construct mul- 
tiple black hole initial data with linear momentum and 
spin (see also [26, 27]). As it turns out, even in the 
general case of orbiting black holes one can analytically 
separate the coordinate singularity at the puncture and 
work on R 3 , both for initial data and evolution [8, 11]. 

Here we concentrate on evolutions of Schwarzschild in 
the puncture framework. The 3-mctric of Schwarzschild 
in isotropic coordinates at t = is 

Wds 2 = {l + ^-) 4 (dx 2 +x 2 (m 2 ), (81) 
2x 

where the radial grid coordinate x coincides with the 
isotropic Schwarzschild coordinate R as given in (39), 
but differs at later times depending on the shift and the 
boundary conditions used during evolution. The coordi- 
nate singularity at the puncture, x = 0, corresponds to 
the inner asymptotically flat infinity of the black hole, 
i.e. the puncture corresponds to the left-hand i° in the 
Carter-Penrose diagrams of Subsec. VI B. 

For evolutions, one can define conformally rescaled 
metric and extrinsic curvature variables, 

9ij = °*- 4 G 13 , k %] = °*- 4 A^-, (82) 

where = 1 + |^ is the time independent conformal 
factor of Schwarzschild in isotropic coordinates. As ar- 
gued in [8] and discussed further in [11], as a consequence 
of this rescaling there are no singularities apart from the 
coordinate singularity in ^, in particular the rescaled 
metric and extrinsic curvature are regular and do not 
evolve at the puncture. What concerns us here is the 
regularity of the maximal slicing condition (12). 

Using the conformal metric gij, one has to solve equa- 
tions of the type 

A 9 a- 0(-)d k a = 0. (83) 

x 

As it turns out, standard numerical methods to solve 
this elliptic equation can find a regular solution for a on 



R 3 , which has the feature that its gradient vanishes suffi- 
ciently rapidly at the puncture such that 0(^)di~a is zero 
at the puncture, and the lapse collapses to zero near the 
puncture. This was noted as an experimental fact in [7] 
and analyzed in [8, 11], but there was no rigorous proof 
that a regular solution to (83) exists, and while numeri- 
cally quite robust, it was not clear that the numerics was 
able to reliably determine this lapse. 

The proposal for the "zgp" lapse is to impose the van- 
ishing of the gradient as a boundary condition on the 
lapse in coordinates that avoid the coordinate singular- 
ity of isotropic coordinates. 

B. Derivation of the puncture lapse 

In order to derive boundary conditions for the "zgp" 
lapse, note that at infinity and at the puncture the radial 
grid coordinate x, and also the isotropic Schwarzschild 
coordinate R at t = given in (39), is related to the 
Schwarzschild radius through r = x for {r^oo,.T^oo} 
and r oc - for {r — > oo, x — > 0}. As before we demand 
unit lapse at infinity, 

lim a+ = lim a+ = 1 Vt Z3p , (84) 

x — >oo r — * oo ^ 

such that the slices are labeled by proper time at infinity. 
At the puncture, we impose the condition of vanishing 
gradient by 

lim d x a~ oc lim r 2 d r a~ = Vr zgp (85) 

x — >0 r — >oc y ^ 

making use of d x oc r 2 d r in this limit. 

It turns out to be convenient to derive the "zgp" lapse 
by making the ansatz [54] 

at gp (C,r)=<S>(C)-at ven (C,r) + *(C)-af dd (C,r) (86) 

for a time-dependent linear combination of the odd and 
the even lapse. This form is possible since the maximal- 
ity condition (12) is linear in the lapse. Furthermore, 
with a dd and a even two linearly independent solutions 
are known as discussed for Schwarzschild in Appendix A 
of [4] and as is true also for Reissner-Nordstrom with 
< |Q| < M. This ansatz reveals useful insights in the 
properties of odd, even, and "zgp" lapse, leading e.g. im- 
mediately to the conjecture of Subsec. VD regarding the 
lapse at the right-hand event horizon at late times. Also, 
the smoothness of the lapse (92) and of the height func- 
tion (97) follow trivially from the smoothness shown for 
odd and even boundary conditions in [4]. Equivalcntly, 
the puncture lapse can of course be derived starting with 
the general solution for the lapse, (27), by applying di- 
rectly the desired boundary conditions (84) and (85). 

Care has to be taken in demanding these boundary 
conditions analytically since condition (84) is imposed 
on the right-hand side of the throat (denoted by a super- 
script "+" ) while (85) on the left-hand side (denoted by 
a "— "). This, however, will be taken into account in the 
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following as the odd lapse is taken in an antisymmetric 
manner across the throat using the appropriate plus and 
minus sign in (29) for the regions extending to infinity 
and puncture, respectively. The even lapse is symmetric 
with respect to the throat by construction. 

Since both of these two lapses are already one at infin- 
ity, the first boundary condition (84), 



$(C) + *(C) = 1 VC, 



(87) 



simply yields that the sum of the dimensionless multi- 
plicator functions <& and \& is one for all times, so one 
can substitute *(C) = 1 - $(C) VC. This (trivial) outer 
boundary condition already allows to formulate a conjec- 
ture regarding the value of the lapse at the right-hand 
event horizon at late times as stated in Subsec. VD. 

Continuing in the derivation of the puncture lapse, the 
gradients of the odd and even lapse at the puncture can 
be obtained by calculating 



lim r 
lim r 



• r2 2 d r a odd 



-M, 



M 



K c (°c)- 



Remembering that the integral Kc{o6) carries units of 
jj and diverges for C ~ * Cn, n , it turns out that for late 
times the gradients at the puncture for a ol id and ot even 
just differ by their sign. Therefore, in order to produce a 
vanishing gradient, at late times the average of the two 
lapse functions has to be taken, i.e. 



lim $(C) = lim *(C) = -. 



(89) 



For arbitrary C, < C < Cu m , we use (88) with the 
second boundary condition (85) to obtain 



*(C) • (M 



K c (oo) 



)-(l-$(C))-M = VC. (90) 



Hence from (87) and (90) the multiplicator functions 
<f> and iff are found to be given by 



$(C) = 1 - <J(C) 



1 K c (oo) 

2 ^ c (oo) + i' 



(91) 



where the divergence of Kq{oo) in the limit C — * Cu m 
gives the average of the odd and even lapse for late times 
as formulated already in (89). 

The "zgp" or puncture lapse now can be written as 



«± ir r \ - ^Pc{r) dC dt zgp 



(92) 



$(C)-at ven (C,r)+V(C)-a± d (C,r) 
1 1 ( 1 



making use of the ansatz (86). Similar to the even case, 
the first line is obtained by passing to Schwarzschild co- 
ordinates, which leads to (45) and (46) with the subfix 
"zgp" instead of "even". 

We show the "zgp" lapse profile in these coordinates 
at times at infinity T zgp = {0, 1M, 2M, AM, 8M} when we 
discuss the numerical evolutions, see Fig. 6. 

The shift and the radial component of the 3-mctric are 
again given by 



and 



Pzgp(C, r) 



lzgp{C, r) 



Qz S p(C,r)7 2ffp (C,r) 



C 



f(r) 



Pc(r) 



(93) 



(94) 



In order to fix the T zsp -dependence of C one may from 
(92) by analogy with (64) conjecture that 



~ \ K c {r)±\Kc{™) + ^ 



(95) 



and hence 
dC'' 



By going back to the differential equation for the max- 
imal lapse, (26), one can readily verify that (95) and 
hence (96) are correct and therefore the time-dependence 
is consistently fixed by the height function 



tf gp (C,r)=Hc(r)±(H c (oc)-^) 



(97) 



using (60) and (62). Furthermore, time is measured at 
right-hand spatial infinity by 



(C)= lim t+ gp (C,r)=2H c (oo) 



c_ 

M' 



Note, however, that the time measured at the puncture 
approaches the finite value 



lim lim t (C,r) 



c 



it in 



M 



(99) 



We are not aware of a physical meaning of this finite time 
at the puncture. 

Finally, we want to point out that using our re- 
sult for the "zgp" lapse, (92), together with r cx - for 
{r — > co, x — > 0}, one can readily verify that the first and 
second derivative of the lapse at the puncture both van- 
ish. So the "zgp" lapse is at least C 2 there. However, 
when calculating higher derivatives of the lapse, terms 
of higher order in the relationship between r and x can 
be of importance, and these terms depend on the chosen 
shift. For a vanishing shift the transformation r(C(r), x) 
will be derived in Subsec. VII C. 
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C. Initial lapse profile 

For the initial lapse profile with (78) the multiplicator 
function $ as in (91) for C = is found to be 



found to approach the value a\ r+ = Cl Jr asymptoti- 
cally. In particular, for the Schwarzschild spacetime with 
°C lim = |V3M 2 
be °a\ r+ = 

Without attempting a complete proof, we can argue as 
follows. Since for fixed time, i.e. for fixed slice label C, 
the lapse arises from the maximality condition (12), i.e. 
a second order linear ODE with smooth coefficients, the 
theory of ODE states that the set of solutions forms a 
2-dimcnsional real linear space, c.f. Theorem I in §15 of 
[33]. As a consequence, every lapse using the superposi- 
tion principle can be constructed as a (time-dependent) 
linear combination of the odd and even lapse, which are 
evidently linearly independent for < \Q\ < M. From 
(29) and (64) evaluated for C — > Cu m , one can see that 
odd and even lapse at the right-hand event horizon in 
this limit arc given by 



$(C = 0) 



M 



M 2 + Q 2 



(100) 



In particular, for Schwarzschild with °$(C = 0) = 1 one 
starts with the even lapse only, whereas for Reissner- 
Nordstrom with increasing charge the balance is shifted 
more and more to the odd lapse to produce a "pre- 
collapsed" lapse profile. 

So the initial puncture lapse profile in terms of the 
Schwarzschild radial coordinate r using (92) is found 
to be given on the right-hand and left-hand side of the 
throat as 



a zgp (C = 0,r) = (M2 + Q2)r ■ (101) 

Applying the transformation (38) one can derive this pro- 
file also for the Schwarzschild isotropic coordinate R, 



at p (C = 0,R) = l 



2Q 2 



AR + M 2 - Q 2 



M 2 + Q 2 4R 2 + AMR + M 2 ~Q 2 ' 

(102) 

Note that for the puncture lapse the isometry (42) is no 
longer present. 



: sj. 

and °r + = 2M this value turns out to 
« 0.3248. 



lim OL^ ven | r+ = 



(103) 



Furthermore, any other lapse constructed by a linear 
combination with coefficients adding up to one in or- 
der to have unit lapse at infinity has the same limiting 
value. Finally, we argue that as the slices approach with 
rc — > 7'Ci im > a limiting surface and hence avoid the 
singularity [5], the limit C — ► Cu m considered so far car- 
ries over to the limit of infinite time, r — * oo. Although 
reasonable, it is not yet clear to us whether the last step 
holds without exception. 



D. Late time limit for the lapse at the right-hand 
event horizon 



VI. COMPARISON OF THE ANALYTICAL 
ODD, EVEN, AND ZGP MAXIMAL SLICES 



Based on the ansatz (86) we want to state a conjecture 
for the lapse at the right-hand event horizon at late times 
that is of some importance for black hole excision. Al- 
though one can use the apparent horizon as the location 
for the excision boundary, finding the apparent horizon 
is usually time consuming, and a simpler prescription for 
some surface that coincides or is located inside the event 
horizon can serve the same purpose for excision. For ex- 
ample, one sometimes can define the excision surface by a 
surface of constant lapse. In numerical simulations with 
maximal slicing, a lapse between 0.3 [28] and 0.34 [29] 
turned out to mark the location of the right-hand event 
horizon. A value in this range can also be inferred from 
some graphs in the literature showing the lapse as a func- 
tion of time, c.f. Fig. 3(a) in [30], Fig. 6 in [31] or Fig. 13 
in [32]. 

For very general boundary conditions we can state the 
following conjecture: 

For maximal slicing of the non-extremal extended 
Reissner-Nordstrom spacetime with the throat of the 
Einstein-Rosen bridge coinciding initially with the 
event horizon and with boundary conditions imposed 
on the lapse other than the odd ones, the lapse 
at late times at the right-hand event horizon is 



A. Embedding diagrams 

Embedding diagrams are useful tools for the visualiza- 
tion of geometric properties of slicings. Following Ap- 
pendix F of [34] and Sec. IV of [35] in suppressing one of 
the angular coordinates, by embedding a 2-dimcnsional 
Ricmannian hypersurface into flat 3D space the intrinsic 
geometry of this surface is preserved. The embedding of 
the Reissner-Nordstrom {r, cf>} line element, i.e. the equa- 
torial hypersurface 6 = ^ corresponding to fixed time at 
infinity t, in Euclidean space written in cylindrical coor- 
dinates {z, r = y 1 x 2 + y 2 , </>} is obtained by demanding 



< 3 W 



dz 2 + dr 2 + r 2 dch 2 



7(r) 



f(r) 



di 



2 + r 2 d(b 2 



pc{r) 



dr 2 +r 2 d6 2 . 



(104) 



For fixed C(r) this results in the differential equation 



dz 
dr 



C— const 



Pc(r) 



(105) 
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for the embedding function z(C, r). For < \Q\ < M and 
f > fc its solution is found by integration, 



z[M] 



z(C,r)=± f 

J rr 



pc(y) 



1 dy, 



(106) 



where without loss of generality symmetry with respect 
to the z = plane has been chosen. In particular for the 
initial Schwarzschild slice, i.e. for Q = C = 0, the solu- 
tion [34] is found to be 



>z(C = 0, r) = 2V2y/(r - 2M)M. 



(107) 



The embedding is a parabola of revolution referred to 
as Einstein- Rosen bridge [17]. Here in both "universes" 
the 3-geometry becomes Euclidean far from the throat, 
which is the minimal 2-sphere of the bridge. 

Note, however, that Einstein's field equations fix only 
the local geometry of the spacetime but not its topology. 
A geometrically identical - but topologically different - 
embedding could be found for example by identifying 
both sheets in order to obtain the "throat of a wormhole" 
in the sense of Misner and Wheeler [23] . Here in one flat 
space two distant regions are connected in the limit when 
the separation of the wormhole mouths is large compared 
to the circumference of the throat, c.f. Fig. 31.5 in [36]. 

For the general < \Q\ < M, < C < Cu m case we 
evaluated the integral (106) numerically using for the 
relationship C(r) equations (52) and (98) for the even 
and "zgp" boundary conditions, respectively. The em- 
beddings of these two cases therefore only differ by their 
time labeling and are less informative than their corre- 
sponding spacetime diagrams shown in Subscc. VI B. In 
particular, information about any underlying symmetry 
of the slices - like the (anti-)symmetry with respect to 
the throat for (odd) even boundary conditions and the 
more complicated behaviour for "zgp" conditions - has 
been lost. The results are shown for Schwarzschild in 
Fig. 2, which should be compared for the even boundary 
conditions with Fig. 4.18 in [34] or Fig. 2 in [3]. Re- 
flecting the curves at the z = plane and rotating them 
around the z-axis, a set of throats is obtained starting 
with a parabola of revolution with the minimal 2-spherc 
initially at the event horizon °r + = 2M and degenerat- 
ing at late times to an infinitely long cylinder with radius 
r = °rc lim = as the slices approach the maximal hy- 
persurface with this radius asymptotically. The plots for 
charge in the range < \Q\ < M, which are not shown 
here, are qualitatively very similar. 



B. Kruskal-Szekeres and Carter-Penrose diagrams 

Spacetime diagrams such as the Kruskal-Szekeres [37] 
and Carter-Penrose [38, 39] diagrams provide a conve- 
nient way of visualizing spacetime properties like causal- 
ity, see e.g. [36, 40]. In this paper we want to limit our- 
selves to the discussion of the Schwarzschild spacetime, 
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FIG. 2: In the embedding diagrams for maximally slicing the 
Schwarzschild metric using odd, even, and "zgp" boundary 
conditions, the Einstein-Rosen bridges are obtained by reflect- 
ing the curves corresponding to r = {0, 1M, 2M, 4M, 8M} 
displayed from the bottom up at z = and rotating them 
around the z-axis. They start with the throat at °r+ = 2M 
and degenerate for r — > oo in the even and "zgp" case - which 
just differ by their time labeling - to an infinitely long cylinder 
"rc Hm = 2M 



the results for the non-extremal Reissner-Nordstrom 
spacetime can be found in [21]. 

For the Schwarzschild metric (32), radial null geodesies 



motivate for r £ 



2M, oo [ the introduction of the 



Regge- Wheeler radial coordinate 
dr 



1 — Ml 



2Mln 



r-2M 



2M 



, e] — oo, oo [. 

(108) 

As G?(t±r*) = on radial null geodesies, ingoing 
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Eddington-Finkclstcin coordinates {v,r, 9, <j>} with 

v = t + r*, ug]-co,oo[ (109) 

leading to ds 2 = -(1 - 2M) cfo 2 + 2 cfodr + r 2 c/X} 2 can 
be defined as in [41]. Similarly, outgoing Eddington- 
Finkelstein coordinates {u, r, 9, </>} given by 

u = t — r*, u e] — co, oo [ (HO) 

lead to ds 2 = -(1 - ^M) du 2 - 2 d-udr + r 2 dQ 2 . These 
line elements can now be analytically extended to r > 0, 
and the corresponding Eddington-Finkelstein diagrams 
are shown e.g. in [40]. In these coordinates the metric 
components are non-singular at °r + = 2M. However, it 
should be mentioned that the regions r < 2M in ingo- 
ing and outgoing Eddington-Finkclstcin coordinates do 
not coincide, since one can readily verify that for all 
timelike or null worldlincs from ds 2 < it follows that 
2dvdr < 0, while 2dudr > 0. Therefore, the nature of 
the singularity at r = is different: An analysis in ingo- 
ing Eddington-Finkelstein coordinates shows that no sig- 
nal can escape back to infinity once it passed the event 
horizon. The time reverse of a black hole, a white hole, 
is found by performing a similar analysis in outgoing 
Eddington-Finkelstein coordinates . 

With the exterior region covered by both in- and out- 
going Eddington-Finkelstein coordinates, one can write 
the Schwarzschild metric as 

2M 

ds 2 = -(1 )dudv + r 2 dn 2 

r 

32 M 3 

= e - ™- dUdV + r 2 dVL 2 (111) 

r 

by introducing for r > 2M the Kruskal-Szekeres coordi- 
nates [37] 

U = -c~™r < o and V = +C+™" > 0, (112) 

which are analytically extendable to U > and V < 0. 
The Schwarzschild radial coordinate r is now given im- 
plicitly as a function of U and V by 

U ■ V = -e* = (l - J_) e 3ir. (113) 
v 2M J v ' 

Note that the singularity at r = corresponds to 
U ■ V = 1 and the event horizon at °r + = 2M to 
U ■ V = 0, i.e. to either U = or V — 0. It is convenient 
to introduce new time and space coordinates by 

T=^(V + U) and X = -(V-U) (114) 

in order to produce Kruskal-Szekeres diagrams. Note 
that with t expressed in terms of T and X by 

1 = 4M { "ctennf in the { II, IV ' (H5) 

the lines of constant t correspond to straight lines 
through the origin, whereas curves of constant r are hy- 
perbolae with asymptotes T = ±X. 



In Fig. 3, we show Kruskal-Szekeres diagrams for odd, 
even, and "zgp" lapse. (Note the difference to the hand- 
drawn Fig. 1 in [11].) Since the in- and outgoing ra- 
dial null geodesies are given by lines of constant U and 
V, the light cones are unit light cones appearing at 45 
degree. One can infer from the Kruskal-Szekeres dia- 
gram that the Schwarzschild geometry consists of four 
regions separated by the (dot-dashed plotted) event hori- 
zon at °r + = 2M, namely the regions I and I', two iden- 
tical but distinct asymptotically flat universes, where 
the isometry °R < — > corresponds to the mapping 
(T,X) < — > (-T,-X), and the regions II and IF, two 
identical but time-reversed regions in which physical sin- 
gularities (a black hole and a white hole) are present. 

As discussed in Sec. Ill, the maximal slices with 
odd boundary conditions correspond to surfaces of con- 
stant Schwarzschild time coordinate t in the outer 
regions I and F since Jj generates a Killing vec- 
tor field there [4, 42]. In particular, in the upper 
plot of Fig. 3 the solid black straight lines in re- 
gions I and F are (starting with the horizontal time- 
symmetric t = slice) the spacelike maximal hypersur- 
faces t = {0, IM,2M,4M,8M}, whereas the gray time- 
like lines correspond to r = {2.5Af, 3M, 3.5M, AM}. In 
regions II and IF the spacclikc limiting maximal slice 
r = °rc lim = ^rf- is plotted as a dashed line, the singu- 
larity r = by a zigzag line. 

In the even case, plotting the height function 
),r) as in (48) for the previously stated 
times at infinity, an observer moves forward in 
Schwarzschild time t in region I and backward in region I' 
in a symmetric manner in order to reach the limiting 
maximal slice, c.f. the plot in the middle of Fig. 3, and 
e.g. Fig. 1 in [3] or Fig. 3 in both [6] and [18]. However, 
in terms of the new time coordinate T even , time runs for- 
ward (due to the always positive lapse) equally fast (due 
to the symmetry) in both regions so that the asymptotic 
values for t at puncture and infinity differ by the amount 
2T even . Due to the symmetry with respect to the T-axis, 
the throat °rc (plotted as an unfilled box) moving from 
°r + = 2M to °rc Um = remains on this axis and the 
slices penetrate into the left-hand and right-hand event 
horizon (denoted by downward and upward pointing tri- 
angles) symmetrically. 

Demanding for "zgp" boundary conditions neither 
symmetry nor antisymmetry, the slices are "lop-sided" 
as can be seen by plotting t Z gp{C(T zgp ), r), (97), in 
the graph on the bottom of Fig. 3. In particular, 
one can observe that as a function of time the throat 
°rc moves "to the right" along the shown curve corre- 
sponding to t = H c (oo) — jj. From (97) one can in- 
fer that as in region I at infinity the time measured as 
T zgp (C) = lim^oo t+ sp (C, r) goes to infinity in the limit 
C — ► Cu, n , as stated in (99) in region P at the punc- 
ture the finite time = |a/3M w 1.2990M is found. 
Taking care of the sign as for even boundary condi- 
tions, the maximal hypcrsurfaces approach in region P 
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In Fig. 4, we show Carter-Penrose diagrams corre- 
sponding to Fig. 3. These are obtained by a conformal 
compactification of the form ds 2 = £l 2 (t,x l ) ■ ds 2 which 
- leaving the underlying causal structure unchanged - 
maps points at infinity in the original metric ds 2 to a 
finite affine parameter in the compactified metric ds 2 . 
Sec [40] for further technical conditions on the spacetime 
and on the (not unique) choice of Q to guarantee that 
this construction will work. As shown in this reference, 
on the conformal boundary one can identify past and fu- 
ture null infinity past and future timelike infinity 
and spacelike infinity i , i.e. the beginning and end of 
null, timelike and spacelike geodesies. 



odd 




iX 



FIG. 3: The Kruskal-Szekeres diagrams for the extended 
Schwarzschild spacetime show the maximal slices for time 
at infinity r = {0, 1M, 2M, AM, 8M} for the odd, even, and 
"zgp" boundary conditions. Note that the puncture corre- 
sponds to the spatial infinity of region I' located to the left of 
the plot. 



FIG. 4: Carter-Penrose diagrams are obtained by compacti- 
fying the Kruskal-Szekeres diagrams of Fig. 3 as stated in the 
text. Note the behaviour of the slices at spatial infinity, the 
movement of the throat rc and of the left-hand and right- 
hand event horizon, which are again plotted as unfilled box, 
downward and upward pointing triangles, respectively. 
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For the Schwarzschild metric, compactification can be 
carried out [42] by introducing as in (114) new time and 
space coordinates T and X based on 



U = arctan U and V = arctan V, 



(116) 



where U and V are in the range — -| < U, V < ^. 

It is worth pointing out that as radial null geodesies 
are again unit light cones the maximal slices arc obvi- 
ously spacclikc. Furthermore, at the right-hand i°, with 
the lapse being one at infinity for all three boundary con- 
ditions, essentially t = const slices are found, whereas in 
particular for "zgp" boundary conditions at the punc- 

the line t = — ^jf^ is a P _ 



ture, i.e. the left-hand 



proached at late times. One may readily verify that for 
the compactification carried out here these slices arrive 
at i° with a slope given by 
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dX 



tanh- 



4Af 



(117) 



In our case, starting with the horizontal time-symmetric 
hypcrsurface t = the slices for all three boundary con- 
ditions fan out to hug J + asymptotically, while for "zgp" 
boundary conditions at the puncture a finite angle given 
by arctan [tanh^p-] sa 17.4° is found in the limit of late 
times. 



VII. COMPARISON WITH NUMERICAL 
SIMULATIONS 

A. Numerical method 

Maximal slicing has been studied numerically for the 
evolution of a single Schwarzschild black hole already in 
[3] , and was later on implemented and used as a standard 
test case in ID- [6], 2D- [31] and 3D-codes [7, 11] based 
on a Cartesian grid and by using a smooth lattice method 
[32]. 

For comparison with the analytical result for the punc- 
ture lapse, we have performed simulations with the 
Cactus code [43]. We evolved a single Schwarzschild 
puncture with maximal slicing and vanishing shift with 
the Baumgarte-Shapiro-Shibata-Nakamura (BSSN) sys- 
tem of evolution equations, see [11] for details. The max- 
imal slicing lapse is obtained by solving the elliptic equa- 
tion (12) using the multigrid solver BAM_Elliptic [44, 45], 
without imposing any additional boundary condition at 
the puncture. 

In Fig. 5, we show an example for the numerical grid. 
The runs are carried out in 3D on an octant of i? 3 , which 
is possible for spherically symmetric systems. The punc- 
ture is located at the origin, and we choose a staggered 
Cartesian grid with uniform grid spacing such that the 
puncture is not part of the grid but rather is located 
half-way between grid points in all three directions. 



throat 



"left- 
hand" 
EH 



punct ure 



Y[M] 
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"infinity" 

"right- 
hand" 
EH 




X[M] 



FIG. 5: The Cartesian grid used in numerical solutions is 
shown (suppressing the Z-direction) for 32 3 grid points and a 
resolution of AX = AY = AZ = 0.05M, which places the 
outer boundary at X t w 1.6M. The domain is one octant 
with appropriate boundary conditions for spherical symmetry 
imposed at the planes X — 0, Y — 0, and Z — 0. Note that 
the puncture is staggered at the origin. The numerically found 
locations of the throat xc and the left-hand and right-hand 



event horizon. 



and 



are superimposed as dark solid 
8M. 



The outer boundary condition for the lapse is imple- 
mented numerically at the edge of the finite numerical 
grid by employing a Robin condition [46] as described in 
[47]. In order to study errors introduced by the Robin 
boundary condition implemented at the edge of the grid, 
runs have been carried out using a uniform spacing of 
AA" = 0.05M and (multiplied by a courant factor of |) 
time steps of At 23P = jAX = 0.0125M for a series of 
32 3 , 64 3 and 128 3 grid points to place the outer boundary 
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at X ob « 1.6M, 3.2M and 6.4M, respectively. Similarly, 
we checked for second order convergence using different 
numbers of grid points for the outer boundary close in. 
These numerical tests are described in [21]. For the res- 
olutions considered here, there are limitations in the nu- 
merical convergence near the puncture, and it is known 
that higher resolution is required to see clean convergence 
near the puncture, sec [11] for algebraic slicings. For the 
present purpose of comparison with the analytical results, 
the goal was only to establish that these runs are approx- 
imately in the convergent regime. All the results in the 
remainder of this section are for AX = 0.05M, 128 3 grid 
points, and X i> « 6AM. 

We will show that the numerical evolution obtained 
in this way agrees nicely (within a certain error) with 
the analytical result for the puncture lapse. In order to 
carry out this comparison, we have to establish an ex- 
plicit mapping between the coordinates used in the nu- 
merical simulation and our analytical study. 



B. Schwarzschild radial coordinate 

One has to remember that for 3D Cartesian coordi- 
nates {X, Y, Z} the system has no spherical symmetry 
from the point of view of the code. However, looking at 
radial output the radial coordinate x defined as the Eu- 
clidean distance from the puncture is readily obtained. 
Because of the spherical symmetry the expression r 2 dfl 2 
for the angular part of the metric is also expected in nu- 
merical simulations [55]. Hence the Schwarzschild radial 
coordinate can be obtained by calculating 



[T,X). 

(118) 

Fig. 6 shows the odd, even, and "zgp" lapse profiles 
for times r = {0, 1M, 2M, AM, 8M } measured at infinity. 
The solid lines show the analytic results based on the 
evaluation of (29), (64) and (92). For the "zgp" bound- 
ary condition, data points obtained from the largest sim- 
ulation arc plotted as boxes in Fig. 6 and lie virtually on 
top of the analytical results. 

Evaluating the analytic integrals is a non-trivial task 
due to the late time divergence of the involved integrals 
Hc(r) and Kc{r) as discussed in more detail in [21] 
based on [48]. The value of the lapse at the throat °rc is 
marked on these slices by an unfilled box and similarly for 
the left-hand and right-hand event horizon correspond- 
ing to = 2M by downward and upward pointing tri- 
angles. Being motivated geometrically these "markers" 
essentially allow one to sketch for a certain time at infin- 
ity the corresponding lapse profile and are hence suited 
for a study of their late time behaviour which will be 
carried out in a later paper [20]. Note, for example, that 
for the collapsing lapse profiles obtained for even and 
"zgp" boundary conditions the radial coordinate of the 
throat approaches °rc lim = and that the conjecture 



of Subsec. V D regarding the lapse at the right-hand event 
horizon holds. 
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FIG. 6: Shown are the odd, even, and "zgp" lapse 
versus the Schwarzschild radial coordinate r for times 
r = {0, 1M, 2M, AM, 8M} measured at infinity. The solid 
lines show the analytic results, while the numerical results for 
the "zgp" lapse are plotted as solid boxes. The numerical data 
points lie virtually on top of the analytical result. Note that 
in the odd case the time-independent lapse is negative on the 
left-hand side of the throat (denoted by an unfilled box) coin- 
ciding with the event horizon, whereas for the even case the 
collapsing lapse is symmetric in the sense that in particular 
left-hand and right-hand event horizon (shown as downward 
and upward pointing triangles) and also infinity and punc- 
ture coincide. These symmetries are lost for the "zgp" lapse 
with the lapse remaining one at infinity and collapsing at the 
puncture. 
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C. Radial isotropic grid coordinates 

From the analytical point of view the task of finding 
the foliation of the extended Reissner-Nordstrom met- 
ric using maximal slices with "zgp" boundary conditions 
has been solved already in Sec. V in the sense that in a 
well-suited coordinate system, namely the radial gauge, 
(13), the maximal slices arising as level sets T zgp = const 
of the height function t zgp (C(T zgp ),r) , (97), have been 
found. Keeping in mind the 3+1 split, the foliations 
shown previously in Kruskal-Szckcrcs and Carter-Penrose 
diagrams are uniquely determined by the lapse func- 
tion af gp (C(T zgp ) 7 r), (92), telling one how to progress 
from slice E r to the consecutive one £ r +,5 T . Since for 
puncture data in Cactus neither the Schwarzschild co- 
ordinates, {t, r, 9, (f)}, nor the previously introduced co- 
ordinates {r, r, 9, (j)}, but depending on the shift differ- 
ent "grid coordinates", {r, x, 9, 4>}, are used, it is for a 
comparison with numerical results important to find the 
coordinate transformations relating these gauges. Above 
all we arc now interested in using the freedom of choos- 
ing different coordinates (x l ) on these slices in order to 
construct coordinates such that the 4-mctric coincides 
at all times with output from a numerical evolution of 
black hole puncture data and hence to prove by construc- 
tion the existence of an underlying analytical solution. 
The radial coordinates of this puncture data line element 
we shall call "isotropic grid coordinates" in the case of 
a vanishing shift and "isothermal grid coordinates" for 
a non-trivial shift to emphasize the connection both to 
the isotropic and isothermal gauge and to the numerical 
grid. Unless stated otherwise, here and in the following 
the expressions are valid for arbitrary boundary condi- 
tions imposed on the lapse. However, one should note 
that for odd boundary conditions the time-independent 
isotropic coordinates of Sec. Ill are recovered, which are 
not horizon-penetrating. In this paper we concentrate on 
the case of vanishing shift, B = 0, where B denotes the 
shift for isotropic coordinates to distinguish it from the 
shift (3 in radial gauge. The more general case of a shift 
B ^ is discussed in [21] and subject of a later paper. 

Independent of the coordinate choice (x l ) on the max- 
imal slices, time at infinity is still measured by r. So the 
task is to find a radial coordinate transformation of the 
form r = r(C(r), a;) to a line element with zero shift B, 



ds 2 



1 



) dr 2 + 2/3 drdr + 7 dr 2 + r 2 dfl 2 



:= -A 2 dr 2 + Gdx 2 



dn 2 



(119) 



Here, = nas been used, and the lapse A and the 

3-mctric Gij with its radial component G as functions of 
t and x have been introduced and are given by 



A(t, x) = u(t, r(r, x)), 
G(t,x) = j{T,r(r,x)) ( ^- 



(120) 



To allow for a direct comparison with numerical puncture 
evolutions as produced by Cactus, the 3-metric is written 
in the form 



{3) ds 2 = Gdx 2 +r 2 dtl 2 



* 4 



qdx 2 -\ r dn 2 



(121) 



obtained by rescaling the 3-mctric with the conformal 
factor 



l + Yx ] 
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Ax 2 



(122) 



Here the coordinate x of puncture evolutions is given 
in terms of the 3D Cartesian coordinates {X, Y, Z} by 
x = \J X 2 + Y 2 + Z 2 , which measures the Euclidean dis- 
tance to the puncture at the origin, and 



9ij(r,x) 



Gi 3 (r, x) 
* 4 (x) 



with g(r,x) 



G{t, x) 
f 4 (x) 



(123) 



has been introduced. 

Again it is worth pointing out that for the initial 
slice the radial grid coordinate x is identical to the 
Schwarzschild isotropic coordinate R denned in (38) and 
hence related to the Schwarzschild radial coordinate r by 



r(r = 0,x) = x^> 2 (x) 



(124) 



as discussed in Sec. III. Note, furthermore, that for 
black hole puncture initial data the 3-mctric Gij is im- 
plemented in isotropic coordinates, hence 

G(t = 0, x) = <J 4 (.t) and g(r = Q,x) = l (125) 

has to hold initially. 

Notice that as a time- independent, trivial case the 
odd boundary conditions are included and lead to the 
isotropic lapse of Sec. Ill together with the formulas 
(124) and (125) valid at all times, which yields the 
static Reissner-Nordstrom metric in isotropic coordi- 
nates. However, due to the negative values of the lapse 
in the region between the puncture and the throat, with 
an observer going backward in time, its numerical imple- 
mentation has been found to be unstable in at least one 
example [49]. 

In the following wc use the observation that the punc- 
ture data line clement for zero shift has to have a 3-mctric 
with a time-independent determinant. As remarked in 
Sec. Ill, this is a consequence of maximal slicing, K = 0, 
in combination with a vanishing shift, B = 0. Consider- 
ing the ansatz (121), one can see that due to the time- 
independence of the conformal factor the determinant of 



the rescaled 3-mctric 



also has to be time- 



independent and hence to be a function of x only. The 
latter can be determined on the initial slice to be given 
by det {g m n} = x 4 sin 2 9 using the initial conditions (124) 
and (125). So it follows that the rescaled radial part of 
the 3-metric can be written as 



j(C(t),x) 



x^ 8 (x) 
r A {C{ T ),x) 



(126) 



18 



It is from G = \t 4 g ; 
for fixed slice label C the ODE 



— = 7 (f§) 2 now trivial to infer 



dr 
dx 



^W(z), 



(127) 



which relates the Schwarzschild radial coordinate r to the 
radial coordinate x for puncture evolutions. This can be 
integrated as 



Uy 



y/pc{y) 



± 



y 2 * 6 (y) dy 



(128) 



using the throat as lower integration limit. Here 
xc denotes the location of the throat on the grid 
as a function of time and the "+" and "— " sign 
applies for the right-hand or left-hand side of the 
throat, respectively. One may readily check that for 
C = the ODE (127) coincides with the ODE (37) 
found for Schwarzschild isotropic coordinates leading 
in that case to x = R with the throat at rc=o = r + 
and xc=o = x + = R + = | y 'M 2 — Q 2 . For later times, 
though, xc has to be found from (128) by demanding 
that the coordinate transformation r = r(C(r), x) is con- 
sistant with the requirement of a vanishing shift. Sec [20] 
for details and for a discussion of the late time behaviour 
of xc in the case of even and "zgp" boundary conditions. 

Given these coordinates, we can finally compare the 
analytically given 4-metric with Cactus output for max- 
imal slicing of a Schwarzschild black hole with "zgp" 
boundary condition and vanishing shift. The time- 
dependent profiles, namely the coordinate transforma- 
tion relating r and x, the "outward moving shoulder" 
obtained for the collapse of the lapse and the growing 
peak in the radial component of the 3-metric, are plotted 
in Fig. 7, 8, and 9, respectively. Similar numerical results 
can also be found in the literature, see e.g. Fig. 3.16 in 
[50] or Figs. 2 and 6 in [31] and Fig. 2 in [7], where in- 
stead of x in the latter references usually a logarithmic 
radial coordinate has been used. 

For Figs. 7 to 9 the reader should bear in mind that 
for odd boundary conditions the time-independent coor- 
dinate transformation relating Schwarzschild radial and 
isotropic coordinates is given by (38), the isotropic lapse 
by (40) and the rescaled radial part of the 3-metric is 
simply g = = 1. Although both the evaluation of the 
underlying analytical expressions (as described in [21]) on 
the one hand and the numerical simulations (with errors 
introduced by finite differencing and the outer boundary) 
on the other hand contain numerical errors, the Cactus 
data points virtually lay on top of the analytic curves. 

We want to emphasize that for arbitrary boundary con- 
ditions the coordinate transformation r = r(C, x), (128), 
holds, with differences in the time-dependence of the 
slice label C(r) and the (in general time-dependent) lo- 
cation of the throat xc- For this - an exception be- 
ing the time-independent odd case discussed previously - 
e.g. the unbounded growth of the radial metric function 



boundary conditions, is not only a feature of the puncture 
lapse. Even worse, in the context of zero shift this slice 
stretching effect is a fundamental property of maximal 
slicings when avoiding singularities as will be discussed 
further in [20]. 

The development of the peak in the metric shown in 
Fig. 9 is the combined effect of local observers falling 
into the black hole and the collapse of the lapse function 
shown in Fig. 8. The coordinate points at smaller val- 
ues of the isotropic coordinate have larger infall speeds, 
causing the radial metric component g to increase to- 
wards smaller x. However, at the same time, there is a 
competing effect due to the use of the singularity avoid- 
ing slicing. The motion of the grid points close to the 
puncture at x = is frozen as the lapse collapses there. 
Well inside, the latter effect dominates, and g cannot in- 
crease in time. This causes the radial metric function to 
develop a peak at a place slightly inside the right-hand 
event horizon. At this location, the difference in the in- 
falling speed of the grid points is large, but the lapse has 
not completely collapsed. 

Finally, we want to mention that if one demands a 
time- independent 3-metric Gy, by looking at the angu- 
lar part of both dj and jij given by r 2 d£l 2 it is ob- 
vious that the coordinate transformation relating r and 
x has to be time-independent and hence to be given by 
(124). Therefore this requirement inevitable leads to the 
isotropic lapse. 
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FIG. 7: The Schwarzschild radial coordinate r as given by 
(128) in the context of "zgp" boundary conditions is shown 
as a function of the isotropic radial coordinate x for maximal 
slices with times at infinity r zgp = {0, 1M, 2M, 4M, 8M}. The 
slices are horizon-penetrating and approach the limiting slice 
r = °rc Hm = ^jf asymptotically. Data points from a Cactus 
run have been plotted as filled boxes, which agree well with 
analytic results. The throat (located as minimum of the radial 
coordinate and depictured by a box) is found to move to the 
right-hand side of its initial location °a;+ = 4r, whereas the 
left-hand and right-hand event horizon (both now correspond- 
ing to r = °r+ = 2M and denoted by downward and upward 
pointing triangles) move to the left and right, respectively. 



g(C,x) 



r 4 (C,x) 



^y), (126), as shown in Fig. 9 for "zgp" 
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FIG. 8: The analytically derived collapse of the lapse is shown 
for "zgp" boundary conditions as solid lines for the slice again 
corresponding to r zgp = {0, 1M, 2M, 4M, 8M}. These curves 
have been obtained by evaluating a^ gp (C, r), (92), in the con- 
text of the coordinate transformation (128) shown in Fig. 7. 
Cactus results from the best resolved run are shown as boxes 
lying almost on top of the analytically found collapsing lapse 
profile. An "outward moving shoulder" is obtained as the 
lapse decreases exponentially to zero at the puncture, at the 
left-hand event horizon, and at the throat, but approaches the 

finite value ^y" = j^V3 « 0.3248 at the right-hand event 

horizon. 
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FIG. 9: For the radial component g of the conformally 
rescaled 3-metric gtj for the "zgp" Schwarzschild case in the 
manner of the previous figures, a comparison of the analytical 
result g(C,x) = ~a7^A > (126), and Cactus output is made. 
Note that whereas g = 1 holds for the odd case, for other 
boundary conditions similar slice stretching effects are present 
in the radial metric component as will be discussed further in 
[20]. 



This observation is relevant for the construction of ini- 
tial data for black hole punctures on a maximal slice. 
Such data can satisfy some of the necessary conditions 
for the existence of helical Killing vectors [51], which 
is appropriate when looking for binary black holes in a 
quasi-equilibrium orbit. As initial condition on the gauge 
one may want to impose not only the condition that the 
lapse is a maximal slicing lapse, but furthermore that it 
is everywhere positive. But as we have just argued, for 
Schwarzschild it is only the isotropic lapse for which the 
metric Gij is time-independent, and conversely we con- 
clude that for positive lapse the metric will not be time- 
independent, which conflicts with the quasi-equilibrium 
condition. This was first noted in [21, 51] and indepen- 
dently in [52]. 



VIII. CONCLUSION AND OUTLOOK 

We have introduced the "zero gradient at the 
puncture" boundary condition for maximal slicing 
of the Schwarzschild and Rcissncr-Nordstrom space- 
times. Comparing analytical and numerical results for 
Schwarzschild, this boundary condition leads to convinc- 
ing agreement with the numerically computed evolutions 
for maximal slicing and vanishing shift. As an appli- 
cation of the analytical form for the "zgp" or puncture 
lapse, we can derive a late time limit for the lapse at the 
right-hand event horizon, which is consistent with previ- 
ous numerical estimates. 

Let us mention several directions for further investiga- 
tions. First of all, in [20] we will give a detailed late time 
analysis for the maximal slicing of Reissner-Nordstrom 
puncture data, extending the results for Schwarzschild 
obtained in [4]. We will discuss there the late time be- 
haviour of the 4-metric at the puncture, the throat and 
the left-hand and right-hand event horizon. Note that all 
of our numerical results refer to the Schwarzschild space- 
time. Future numerical work could include an electric 
charge in the range <| Q \< M. 

Maximal slicing is a special case of constant mean cur- 
vature slicings, see for example [19]. It would be inter- 
esting to see whether constant mean curvature slicing is 
amenable to a similar analysis as presented here, and 
whether there is a promising avenue for numerical simu- 
lations. 

As mentioned in the introduction, maximal slicing has 
in some cases been replaced with algebraic slicings of the 
"1+log" type, see [11] for recent work on puncture evo- 
lutions. It may be possible to gain analytic insight into 
1+log slicing, which in fact in some regards mimics max- 
imal slicing. For example, one would like to understand 
the singularity avoidance properties of 1+log slicing bet- 
ter. A first attempt to study this problem has met with 
technical difficulties, but further work is certainly war- 
ranted. 
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Of particular interest from the point of view of numer- 
ical relativity is an analytical study of elliptic shift con- 
ditions in Schwarzschild, say the minimal distortion shift 
condition. A maximal slicing lapse together with a mini- 
mal distortion shift have been suggested as a natural co- 
ordinate choice for numerical simulations [1]. While the 
minimal distortion shift has not been successful for punc- 
ture evolutions, a closely related shift condition called 
Gamma freezing or conformal 3-harmonic has to some 
extent solved the long standing slice stretching problem 
of maximal slicing for single black hole and head-on col- 
lision simulations [11]. In that reference this shift is im- 
plemented as an evolution equation for the shift, but it is 
the elliptic version which is more easily studied analyti- 
cally. A key feature of the Gamma freezing shift is that 
when it is combined with maximal or 1+log slicing an 
approximately time independent metric is obtained for 
the final, static black hole. While in the present paper 
we answered the question whether there exists a maximal 
slicing corresponding to the numerically obtained "zgp" 
lapse, we now can ask the question whether for "zgp" 



maximal slices there exists a shift condition such that 
the evolution of the metric in the resulting coordinates is 
minimized. Some steps in this direction can already be 
found in [21]. 
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